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We performed extensive simulations accompanied by a detailed study of a two-segment size ran- 
dom sequential model on the line. We followed the kinetics towards the jamming state, but we paid 
particular attention to the characterization of the jamming state structure. In particular, we stud- 
ied the effect of the size ratio on the mean-gap size, the gap-size dispersion, gap-size skewness, and 
gap-size kurtosis at the jamming state. We also analyzed the above quantities for the four possible 
segment-to-segment gap types. We ranged the values of the size ratio from one to twenty. In the 
limit of a size ratio of one, one recovers the classical car-parking problem. We observed that at low 
size ratios the jamming state is constituted by short streaks of small and large segments, while at 
high values of the size ratio the jamming state structure is formed by long streaks of small segments 
separated by a single large segment. This view of the jamming state structure as a function of 
the size ratio is supported by the various measured quantities. The present work can help provide 
insight, for example, on how to minimize the interparticle distance or minimize fluctuations around 
the mean particle-to-particle distance. 
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INTRODUCTION 



Since the 1939 paper by Flory [j, for the deposi- 
tion of dimers on lattices, and the 1958 paper by Renyi 
0, 0, for the deposition of segments on a line, the 
random sequential adsorption (RSA) model has become 
a paradigm for the study of many natural phenomena, 
not only in the traditional area of physical chemistry (re- 
action in polymer chains, chcmisorption, colloids, etc), 
but also in the less traditional areas such as biology, 
ecology, and sociology [Jj Ll 0, S El H3- For instance, 
the use of RSA to determine the limiting coverage on 
surfaces requires an uniform deposition from stabilized, 
diluted suspensions of particles sized from 100 A up to 
a micron by convective flow ^(|- Actual objects in- 
volve proteins and submicron colloidal particles [llj . Re- 
cently, experimental interest has also included deposition 
on p atterned surfaces, prepared by lithographic methods 
12], and sonic theoretical work has been also performed 
|l3l Il4| . The basic dimer deposition model has suffered 
several extensions and generalizations, namely, coopera- 
tive sequential adsorption with adsorption rates depen- 
dent on the local environment inclusion of relax- 
ational mechanisms such as detachment and dif- 

fusion [13, Gil , and multi-layer deposition [ljj, |2(| . Also, 
the RSA model and its extensions have been studied in 

one EHLIMIliHEpJ^^ 

and two dimensions [25|, [2(| I27L l29l I30L 13a , l36l l37| either 
in the continuum and lattice versions. More extended 
accounts can be found on recent reviews by Privman 
01111113 and Evans @. 
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Quite recently, interest in the field has shifted towards 
the competitive deposition of mixtures of segments with 
different sizes on the line, with some controversial results 
|32l 1.3 1| . Apart from issues concerning the competitive 
deposition and consequent adsorption of particles at in- 
terfaces, there is a strong motivation for the study of the 
resulting "patterned" structure, either induced by the 
kinetics of deposition |13( or by more controlled means, 
e.g. , by patterning the surface available for deposition 
|l4|. In this paper, we focus our interest on the study 
and characterization of the inter-particle distance distri- 
bution functions of binary mixtures in one dimension. 
We were, thus, able to perform a more detailed study 
of subtle correlations developed during deposition, by 
measuring not only the time dependence of the coverage, 
but also, more refined quantities, such as the distribution 
functions of the distance between particles. From these 
basic measurements of particle-to-particle distance at the 
jamming state, we studied the size ratio dependence of 
the first four cumulants. More specifically, since the third 
and fourth cumulants are straightforwardly related to the 
skewness and kurtosis, respectively, we actually used the 
latter quantities in order to characterize the gap-size dis- 
tribution functions. We observed non-trivial effects, as 
we varied the size ratio of the segments being deposited, 
and we were able to explain phenomenologically some 
of the qualitative features observed in our simulation. 
Therefore, it is justifiable to put some effort to under- 
stand such a clean study case, where one can learn the 
actual effects leading to the cooperative behavior induced 
by deposition towards the jamming state. Moreover, de- 
spite our model being one-dimensional, our analysis can 
serve as a guide to interpret and/or compare with simi- 
lar results in higher dimensions. The paper is organized 
as follows: we present the model and the particulars of 
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the simulations in Section [n] In Section ITTT1 we present 
our results and discuss them. Finally, we present our 
conclusions in Section llVl 



II. THE MODEL AND THEORY 

We consider the competitive deposition on the line of 
segments of two different sizes, namely short segments, 
which we denote as A-segments, and long ones, which 
we denote as i?-segments, under the condition that they 
must not overlap with each other upon adsorption, and 
therefore, mimic an excluded volume, short-range inter- 
action. The fraction of the line occupied by the adsorbed 
segments defines the coverage. We notice that it is possi- 
ble to re-scale, without loss of generality, the length scale 
of the system so that ^4-segments are of unit length, while 
the size of i?-segments is R. In this respect, one can re- 
gard the ratio of the length of _B-segments relatively to 
the length of ^4-segments as the size ratio, 



R 



Length of a B-segment 
Length of an A-segment 



(1) 



The deposition flux represents the number of incoming 
segments per unit length (in one dimension) and per unit 
time. Let us denote by a the corresponding deposition 
flux of A-segments and by (3 that of i?-segments, there- 
fore, having a total incoming flux of segments, a + (3. 
The probability of having an v4-segment attempting de- 
position on the line during an interval of time dt is 



PA 



a + (3 ' 
while that of a B-segment is 



(2) 



PB 



ft 



a + f3 ' 



(3) 



i.e., pb = 1 — Pa as expected. 



Adsorption on the line can only take place if the in- 
coming particle does not overlap with a previously ad- 
sorbed segment, thus mimicking an excluded volume, 
short-range interaction. To make a more straightforward 
comparison with experimental results, we measure time 
in terms of the number of layers of segments, which at- 
tempted deposition whether these segments actually ad- 
sorb on the substrate. 



We performed a series of Monte Carlo simulations for 
various values of the size ratio and equal fluxes of incom- 
ing segments. To reduce the uncertainties in the various 
quantities, it is more relevant to increase the size of the 
system than the number of samples. Therefore, to char- 
acterize the jamming state, we simulated a system size of 
10 7 units and generated 10 2 samples. Whenever we fol- 
lowed the time dependence, we used a different algorithm 
to simulate, which is computationally more demanding 
both in time and allocated memory. Consequently, the 
simulated system size is smaller, 10 4 , and we let time 
evolve up to 100 units, and obtained 10 2 samples. As 
soon as the total size of the segments, which attempt ad- 
sorption on the line, equals the size of the system, we 
increase time by one unit, regardless of the fact that the 
segments actually adsorb, or not. 



We now present some useful definitions and relations. 
For the sake of simplicity, the definitions are valid for 
the jamming state. However, extension to intermediate 
states at specific times only requires the explicit time 
dependence taken into account and also set the upper 
limit of integration to infinity to account for gap sizes of 
all lengths. The probability distribution of empty space 
is defined as 



P%{x)dx = 



Number of empty intervals of size x within da; 
Total number of empty intervals 

I 



da; 



(4) 



in the limit of dx being an infinitesimal quantity and the 
number of ensembles going to infinity. Thus, P$(x)dx for 
the jamming state has the property 



/ P % (x)dx = 1 . 
Jo 



(5) 



Therefore, the jamming coverage is given by, 



yj = 



1 — / xP®{x)dx = 1— < x >0 
Jo 



(6) 



Discriminating all gaps between pairs of consecutively 
adsorbed segments by A A, AB, BA, and BB and 
defining density distribution functions Pb(x), with b G 
{AA, AB, BA, BB}, one obtains the relation 
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P % (x) = P AA (x) + 2P AB (x) + P BB (x) , 



(7) 



where we exploited the fact P A b(x) = Pba{x) for the 
present random sequential adsorption model. Note from 
equations (JSJ and Q that the P&(x) are not normalized. 
Keeping in mind the above definitions, it is now straight- 
forward to reckon higher moments of the gap-size distri- 
bution functions, defined by 



Jtx"P a (x)dx 

Jo Pa(x)dx 



(8) 



with a 6 {0, AA, AS, BA, BB}. Using equation one 
can relate the moments of the gap distribution functions 
with the corresponding moment of the global distribution 
function given by equation © yielding 



< / > = < x 11 > AA [ P AA (x)dx 
Jo 

+2 < / > AB / PAs(x)dx 
Jo 

;(x)dx , (9) 



< X" >BB 



defined for all values of n = 0, 1, 2, ... We also compute 
the cumulants, k^, of a distribution function defined as, 



lnG a (fc) = 



ml 



(10) 



where G „ (k ) is the so called characteristic function de- 
fined by |2i 



G a (k) =< e «- > a = ^t^t^ • W 
Jo ^a(x)da; 

Therefore, from equations JHJ, H10|l . and 111(1 . one derives 
the first four cumulants as 



TABLE I: The table presents the coverage, 9j, as a function 
of the size ratio, R, while a represents the associated error. 
These results are for a system size of 10 7 and for 10 2 samples. 



R 


0j 


CT 


1.00 


0.7475958 


0.000067 


1.05 


0.7544753 


0.000062 


1.10 


0.7599829 


0.000063 


1.20 


0.7683652 


0.000058 


1.25 


0.7716358 


0.000063 


1.30 


0.7744377 


0.000054 


1.40 


0.7789946 


0.000063 


1.50 


0.7825520 


0.000064 


1.60 


0.7854349 


0.000053 


1.75 


0.7890016 


0.000066 


2.00 


0.7941038 


0.000058 


2.10 


0.7961414 


0.000053 


2.25 


0.7991661 


0.000058 


2.50 


0.8036444 


0.000060 


2.75 


0.8073996 


0.000065 


3.00 


0.8105172 


0.000066 


3.25 


0.8131850 


0.000060 


3.50 


0.8155590 


0.000065 


3.75 


0.8176675 


0.000073 


4.00 


0.8195500 


0.000070 


4.50 


0.8227540 


0.000063 


5.00 


0.8253667 


0.000067 


5.50 


0.8275848 


0.000074 


6.00 


0.8294481 


0.000076 


7.00 


0.8324537 


0.000076 


8.00 


0.8347729 


0.000094 


10.00 


0.8381112 


0.000097 


12.00 


0.8403695 


0.00010 


14.00 


0.8420603 


0.00011 


16.00 


0.8433187 


0.00013 


18.00 


0.8443186 


0.00012 


20.00 


0.8451377 


0.00012 



'n 


= < 


X ^> a , 




(12) 


k 2 


= < 


X 2 >a - 


■ <x>l , 


(13) 




= < 


X 3 > a - 


-3 < X 2 > a < X > a 








+ 2 < 


x>l , 


(14) 


ft' 4 


= < 


X 4 >a - 


4 < X 3 > a < X > a 








-3 < 


x 2 > 2 a +12 < x 2 > a < x >l 










-6<x>i , 


(15) 



where n\ is just the mean value and k 2 is the variance. 
The third and fourth cumulants are used in the definition 
of the skewness, 



and kurtosis, 



respectively [36 



S a = 



K a 



a\3/2 



tor 



(16) 



(17) 
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III. RESULTS AND DISCUSSION 

We start by presenting our results with a study of the 
coverage for the particular case of equal fluxes of incom- 
ing segments. In this respect, our study represents the 
particular case of q = 1/2 of reference 34], but with 
a wider range of size ratios. In fact, results presented in 
[S3] are valid for values of R < 2. Our study also includes 
quantities, such as the dispersion, the fraction of empty 
spaces, etc., which are not possible to be computed by 
their method. In Fig.^a), we show the time dependence 
of the coverage, 9(t), for some of the values we simulated 
of the size ratio, R. However, we performed simulations 
for a wider range of values of the size ratio, namely for 
values of 1, 1.05, 1.1, 1.2, 1.25, 1.3, 1.4, 1.5, 1.6, 1.75, 2, 
2.1, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4, 4.5, 5, 5.5, 6, 7, 
8, 10, 12, 14, 16, 18, and 20. Notice that our model for a 
size ratio of one boils down to deposition of segments of 
unit size on the line. The value of 74.75958% ± 0.0067% 
for the jamming coverage is quite close to 74.75979202%, 
first obtained by Renyi and subsequently calculated with 
larger precision by Blaisdell and Solomon 0, H^, H3 • E3 

The jamming coverage increases monotonically as a 
function of the size ratio, as shown in Fig. ^b) and Ta- 
ble[I] This result is not entirely surprising, since for equal 
relative fluxes of incoming segments and for large, asymp- 
totic values of the size ratio, a stretch of the line is either 
fully covered by a large segment or paved by the small, 
unit size, segments with an upper limit of the coverage 
given by the above Renyi value. As both segment sizes 
can attempt deposition with equal probability, we get 
1/2(1 + 0.7476) = 87.38%. However, this upper limit is 
not attained as one increases the value of the size ratio: 
the limiting value of the coverage around 84.5%, obtained 
for a size ratio of 20, remains lower as one can observe in 

Fig.ntb). 

In order to better analyze the structure of the jam- 
ming state, we now proceed to characterize the fraction 
of empty space distributed in terms of the different pairs 
of segments as we vary the size ratio. In Fig. |2Ja), we 
plot the fraction of the substrate left empty at the jam- 
ming state as a function of the size ratio for each pair 
of consecutive segments, namely of type AA, AB, and 
BB. In Fig.[2Ib), we also present the corresponding nor- 
malized population of the gap types. Consequently, in 
the latter case, the sum of the contributions of the four 
gap types adds up to the unit for every value of the size 
ratio. We observe that the fraction of empty space due 
to the AA-gaps increases monotonically with the size ra- 
tio. This can be understood by considering the signifi- 
cant drop in the population of -B-B-gaps at the jamming 
state as the size ratio increases, since these gaps must be 
smaller than unit. The probability of having two consec- 
utive B-segments adsorbed on the line rapidly decreases 
with the size ratio. Not only does the adsorption of both 
large segments imply a clean substrate, for example by 
the absence of smaller segments previously adsorbed on 
the line, but also both large segments must adsorb at a 



-0.5 




0.75 



0.72 1 1 1 ' 1 1 

4 8 12 16 20 

Size Ratio 

FIG. 1: Coverage dependence on the size ratio for equal de- 
positing fluxes of each segment size, a) Plot of the coverage up 
to 100 time units for various values of the size ratio, namely, 
1, 1.1, 2, 4, and 10. b) Plot showing the jamming coverage 
dependence on the size ratio. 



distance smaller than unit to prevent an A-segment to 
fit in. The latter situation becomes less and less proba- 
ble as the size ratio increases, and, to reach the jamming 
state, the leftover empty space between these large seg- 
ments must be hlled with segments of unit size, thus, con- 
tributing to a higher fraction of empty space associated 
with AA-gaps and, consequently, to a higher value of the 
coverage. The slower decrease of the A-B-gap population 
and of the corresponding coverage indicates a jamming 
state constituted by alternating streaks of A-segments 
separated by a single -B-segment. This argument is fur- 
ther substantiated by the strong decay of the fraction of 
empty space associated with the -B-B-gaps for increasing 
values of the size ratio. This is also corroborated by the 
smoother decay of both the fraction of empty space and of 
the population of the A-B-gap type, due to the breakdown 
of the -B-B-gap type as the size ratio increases. Finally, 
it can be seen in Fig.^b) that as the size ratio diverges, 
the relative population of AB- and -B-B-gap types as com- 
pared to the corresponding population of the AA-gaps 
becomes less and less relevant, therefore leaving, in this 
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1 2 3 4 5 0.2 0.4 0.6 0.8 

Size Ratio Size Ratio 



FIG. 2: In both plots, one has the AA-, AB-, and SB-gap 
types represented by circles, squares, and triangles, respec- 
tively, a) Fraction of available empty space, at the jamming 
state, b) Normalized population of gaps at the jamming state. 
In the insets the size ratio varies between 1 and 20. Please, 
refer to the text for further details. 



FIG. 3: Gap-size distribution functions at the jamming state: 
a) For a size ratio of one, the solid curve represents P$(x), 
while the remaining AA-, AB-, and SB-gap types are identi- 
cal, and therefore fall onto a single curve, b) For a size ratio 
of 2.0, one observes the splitting of the various gap types. 
Please, refer to the text for further details. 



limit, long streaks of A-segments. 

In Fig.ljja), we show the distribution functions, P$(x), 
Paa(x), Pab{x), and Pbb(x), at the jamming state, for 
a size ratio of one. Of course, at this size ratio (R = 
1), all gap types are equal as there is a single segment 
size. Therefore, one expects, as shown in Fig. OU^a), to 
observe the collapse of the Paa(x), Pab(x), and Pbb{x), 
onto a single curve. The Pq(x) distribution function also 
satisfies the relation P$(x)/4: = Paa(x) = Pab{x) — 
Pb b (x) in agreement with equation (JZJ) for a size ratio of 
one. In the case of part b) of the same figure, we merely 
changed the value of the size ratio to 2.0, thus breaking 
the symmetries between the various gap types. Once 
again, adding the various gap distribution functions, for 
every value of x, accordingly to equation Q, actually 
reproduces P^{x). As compared to part a), the AA-gap 
type increases its dominance, as pointed out above, while 
the AB-ga,p slightly lowers its contribution. The BB-g&p 
significantly drops its influence, because it becomes less 
probable to have consecutive deposition of B-segments 



with gap lengths smaller than unit. 

Now, wc proceed to the analysis of quantities involv- 
ing cumulants up to the fourth order, more specifically, 
the mean distance, dispersion, skewness, and kurtosis, 
as defined in equations i|12ll7fl . For each of the latter 
quantities, we studied their dependence on the size ratio 
at the jamming coverage limit as shown in Fig. Both 
the AA- and AB-g&j) types have rich non-monotonic be- 
havior, for the mean gap size, dispersion, skewness and 
kurtosis. For example, the minimum values of the gap 
size for the AA- and ^4-B-gaps are not coincidental, as the 
minimum of the AA-gap occurs at a slightly higher value 
of the size ratio. The behavior of the dispersion shows 
the presence of both a minimum and maximum values. 
However, the opposite happens for the maximum value: 
the value of the size ratio at which it occurs is slightly 
lower for the AA-gaj> type. The minimum value occurs at 
the same value of R for these two gap types, within the 
error bounds. The skewness shows a maximum occur- 
ring at the different R values, namely, 1.34 and 1.29 for 
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FIG. 4: Plots involving cumulants up to the fourth order 
of the gap-size distribution functions for each gap type as a 
function of the size ratio. We use the same legend as in Fig. [5] 
a) Average distance between pairs of segments, b) Dispersion 
of the distance between pairs of segments, c) Skewness. d) 
Kurtosis. 



TABLE II: The table summarizes the estimated values of R, 
where the first four cumulants take maximum and minimum 
values, all measured with an error of ±0.005. The correspond- 
ing values of the cumulants are also included. 



gap types 



minimum 
R Value 



maximum 
R Value 



Distance 



Dispersion 



Skewness 



Kurtosis 



AA 
AB 
BB 

AA 
AB 
BB 

AA 
AB 
BB 

AA 
AB 
BB 



1.31 0.2924895 
1.22 0.3291302 

(a) 

1.55 0.2681662 
1.55 0.2788692 
1.01 0.2815730 



1.02 
1.05 



1.34 
1.29 



1.40 
1.36 



0.2825507 
0.2825593 



0.9227950 

0.7253401 


-0.2290402 
-0.6281340 

m 



"The SB-gap is strictly monotonic. 



the AA- and AB-g&p types, respectively. The kurtosis 
reveals a single maximum for these two gap types, once 
again non-coincidental in their size ratio value, with the 
maximum of the size ratio of the AA gap slightly above, 
than the corresponding one found for the AB gap. All 
these values of R are summarized in Table [H] Regard- 
ing the BB-g&p type, one finds a less rich behavior as 
compared with the other gap types. For this gap type 
the mean gap size strictly increases monotonically with 
R, while it decreases strictly monotonically for both the 
skewness and kurtosis. The exception to this behavior 
for the BB gap regards the dispersion, where it shows a 
minimum for R = 1.01 as shown in TablelTTIand Fig.^Jb). 
Finally the dispersion of both the BB- and AB-g&p types 
intersect at R — 6, as shown in the inset of Fig. Hfb). 

If one were to solely deposit large segments, then all 
gap sizes smaller than R would be present. Since we 
are concomitantly depositing unit segments, all gap sizes 
larger than unit must disappear at the jamming state. 
As soon as the typical size of gaps drops below R, the 
large segments stop adsorbing on the line. The smaller, 
unit segments, contrary to larger segments, have to fit 
into all available space, i.e., in all gaps larger than unit. 
For size ratios R < 1.55 the population of BB-g&ps re- 
mains significant, as compared to the AA-g&p one, with 
the consequent formation of B-segments gap sizes close 
to unit. Some of these gaps are really close to unit, but 
strictly larger than unit, in size, and we call these events 
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snug fits. As R increases the number of -BB-gaps de- 
creases, to form A A- and AB-gaps, therefore, increasing 
the probability of snug fit events. The probability that 
for large values of R B£>-gaps form a snug fit also de- 
creases. Consequently, the importance of snug fits de- 
creases for large values of R in agreement with Fig. 
where one observes monotonic behavior of the various 
cumulants in this regime. However, the large number of 
these events for values of R < 1.55 accounts for such rich 
behavior of the cumulants up to the fourth order. For 
example, the existence of minimum values of the AA- 
and A_B-gaps sizes can now be understood as events from 
late stage kinetics close to the jamming state. The same 
argument also applies to the minimum values of the dis- 
persion, since this events tend to lower the uncertainties 
of adsorbed segments. More suble, it is the presence of 
maximum values for the both the AA- and A-B-gaps and 
a minimum of the BB-gaps for values of R < 1.05. Since 
the difference in size of both segments is small, both seg- 
ment sizes compete until coverage values become close 
to the jamming state, thus leading to snug fit events of 
the BB-g&p. This competition leads to lessen the uncer- 
tainty in the -BB-gaps, but it increases the uncertainties 
of the A A- and AB-g&p types. For values of R > 1.55, 
the significant drop in the population of £?_B-gaps can be 
understood of the early onset of a mean gap size smaller 
than R at low coverage values, which effectively blocks 
the adsorption of large segments. Consequently, one ex- 
pects a flatter distribution function of the BB-gap type 
with increasing values of the size ratio, i.e., the net effect 
of increasing the size ratio is to diminish the asymme- 
try of the BB-gap distribution function. The monotonic 
behavior of the BB-gap for large size ratios (R > 1.55) 
stems from the absence of a competing mechanism. Fi- 
nally, the non-monotonic behavior of the remaining gap 
types is due to snug fit events, which the latter tend to 
askew the corresponding gap-size distribution functions 
by favoring small gap sizes. 

The jamming state for large size ratios, follows, there- 
fore, the picture of alternating streaks of ^-segments in- 
terfaced with a single B-segment for large values of the 
size ratio. The presence of such streaks of A-segments 
interfaced with a single B-segment, prevents the AB-g&p 
population to decrease slower than that of the -BB-gap 



one. The snug fits events also tend to favor small gap 
sizes, and this effect makes the distribution function less 
flat and more askewed, thus increasing the value of the 
kurtosis and skewness of both the AA- and AB-gaps. 
(Fig-Hc-d) and Fig.© 



IV. CONCLUSION 

We analyzed in detail the jamming structure of a model 
of random sequential adsorption on the line with two- 
segment sizes depositing with equal fluxes. The struc- 
ture of the jamming state is determined by the popula- 
tion of the various gap types, namely, the AA-, AB-, and 
BB-gap types. For large values of the size ratio, i.e., for 
size ratios greater than 1.55, snug fits become less signifi- 
cant, and the average distance, dispersion, skewness, and 
kurtosis monotonically approach their asymptotic values. 
For values of the size ratio below two, the behavior of the 
above quantities for the AA-, and AB- types are non- 
monotonic due to the presence of snug fit events. The 
jamming state at values of R above 1.55 is characterized 
by streaks of A- interrupted by a single B-segment. At 
sizes ratios smaller than 1.55 a rich, non-monotonic be- 
havior of the above quantities plays develops due to the 
interplay provided by snug fits. From an experimental 
point of view, the study provides to some extend, among 
others, the insight on how to tune up the mean interpar- 
ticle distance or how to minimize the flutuactions of the 
interparticle distance around the mean value. 
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